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Abstract 

In this article we present the structure of the F tests, the variance components and the approxi¬ 
mate degrees of freedom for each of the eight possible mixed models of the strip-split plot design. 
We present an example to illustrate the model and compare it to more traditional settings like 
a three-way factorial design and a split-split plot model. 
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1 Introduction and Method 


There are many opportunities in which a researcher needs to know the behavior of a factor in 
relation to one and/or two additional factors. When this happens it is usually considered a split 
plot design adding subplots into the former plots. It is due, in part, to the great development 
reached by this type of experiment. 

Instead, we study here the strip-split-plot design; i.e., an extension of strip-block designs such 
that each plot on the intersection is subdivided into subplots to insert a third factor. This new factor 
will be more precise on its measurement due to its high number of observations and interactions; 
which is the more important feature of the design. 


We do not claim originalit y on the invention of t his model. On the contrary, iGomez and Gomez 


( 1984l l described it, as well as Zimmermann ( 2004 1 did. They also described the F tests when the 
effects are fixed. Nonetheless, after an intensive search, we could not find on the literature those 
same F tests for the strip-split plot design with mixed effects. Our aim here is to fill that gap. On 
a sequel, we will consider the contrasts for this design and construct their variances and variance 
estimators, again for every case of the mixed effects model. 


To determine the variance components and the ANOVA, we use a method explained bv ISearle et. al. 
( 1992I ). The design is completely randomized so that it makes sense to implement F tests. The 
mathematical model is given by 


Uhijk — m + Rh + Ai eAf^i + Bj + + ABij -\- -|- ACik + BCjk + ABCijk + cj, 

where m is the general mean, is the h-th random block effect {h = 1,... ,r). A* is the i-th. 
horizontal strip effect (i = 1,...,a), Bj is the j-th vertical strip effect {j = 1,... ,b) and Ck is the 
fe-th effect on the intersection of A and B (fe = 1,..., c). So yhijk represents the observation of the 
z-th level of A, the j-th level of B, the k-th level of C on the block h. The errors , eAB^ij 
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and are normally distributed with mean zero and variance cr^, cr^, and cr^, respectively. 
Since the blocks are random, we will assume R iV(0,4). 

The analysis is done according to the scheme on Table [H where df stands for degrees of freedom 
and SS stands for the sum of squares of the respective variation source. 


Table 1: Sums of squares and Degrees of freedom 


Source 

df 

SS 

R 

r — 1 

abcYfh=iiyh... -y...y 

A 

a — 1 

hcrY.tiiv.r..-V..y 

eA 

(r- l)(a- 1) 

bcYft=iYJh=i{yM.. -Vh... -y.L. + y...f 

B 

6-1 


es 

(r- 1)(6- 1) 

- Vh... -y..o. + y...y 

AB 

(a - 1)(6- 1) 

c?'ELi Ei=i (y.ij. -y.i.. -y..j. + y...y 

Gab 

(a - 1)(6- l)(r - 1) 

yhi.. yh.j. y.ij.~^ 

Vh... + y.z.. + -y....)‘^ 

C 

c — 1 

abrJ2k=iiy...k-y....)'^ 

AC 

(a - l)(c- 1) 

brJ2i=lJ2l^l{y.,.k-y.^.. -y...k + y... f 

BC 

(6-l)(c-l) 

(^i'T,’’j=iT,k=iiy..jk -y..j. -y...k + y....)^ 

ABC 

(a - 1)(6- l)(c- 1) 

^1 iy.ijk y.i.k y..jk ~^y...k 

y.ij. ^y.i.. ^y..j. -y....f 

et 

ab{c — l)(r — 1) 

f _ _ j_ 

A^i=l A^j=l 2 -^k=l\yhijk ~y.ijk ~yhij. + 

y.^,? 


2 Sums of squares and covariance matrices as Kronecker products 


Following the algorithm proposed bv iMoser and Sawyer (jl998l ). the sums of squares and the covari¬ 
ance matrices can be computed using Kronecker products. Given the usefulness of this approach 
for computational purposes, we present those products here. Then the sums of squares can be 
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expressed as: 


SSr = Y'[{Ir - r-^Jr) ^ a-^Ja 0 b-^Jb ® c-^Jc\Y, 

SSa = Y'[r-^Jr ® {la - a-^Ja) ® b-^Jb ® c-^Jc]Y, 

SSe^ = Y'[{Ir - r-^Jr) ® {la - Ja) ® Jft ® C"^ 

SSb = Y'[r-^Jr ® a-^Ja 0 {Ib - 6“^ Jfe) ® c-^Jc]Y, 

SSes = Y'[{Ir - r-^Jr) ® ® {Ib - b-^Jb) 0 Jc]l", 

SSab = Y'[r-^Jr ® {la - a-^Ja) ® {Ib - b-^Jb) ® C-^Jc]Y, 

SSe^^ = Y'[{Ir - r-^Jr) ^ {la - a"^ Ja) ® {h ” 5"^ Jfc) ® C"^ J,]!", 
SSc = Y'[r-^Jr 0 a-^Ja 0 6“^ Jft ® (Ic - c"^ Jc)]l", 

SSaC = Y'[r-^Jr 0 (/a - a-^Ja) ® Jfe 0 (/^ - c"^ , 

5Sbc = Y'[r-^Jr 0 a-^ Ja 0 (/^ - b-^Jb) ® (Ic - c-^Jc)]Y, 

SSabc = Y'[r-^Jr 0 {la - a-^Ja) ® (Ife - b-^Jb) 0 (/c - C-^Jc)]Y, 
SSe, = Y'[Ir ® {la - a-^ Ja) ® {h - 6'^ Jfe) ® Ic]Y, 

where Jj is a squared matrix of ones and size j, and /* is an identity matrix of size i. 
The covariance matrix when all effects are fixed is: 

V{Y) = aj^{Ir ® Ja^Jb® Jc) + O-gA {Ir® la® Jb® J c) + 

(^Ib(I^ ®Ja®Ib® Jc) + ® la® Ib® Jc) + 

Crl^{Ir ®Ia®Ib® Ic)- 

The covariance matrix when all effects are random is: 

V(Y) = aR{Ir ® Ja®Jb® Jc) + (^A{Jr ®Ia®Jb® Jc) + 

(^l^{Ir ® la® Jb® Jc) + CrB{Jr ®Ja®Ib® Jc) + 

(^lB{Ir ®Ja®Ib® Jc) + (^AB{Jr ®Ia®Ib® Jc) + 

{I r ® I a ® Ib ® J c) ® ^c{J c ® J a ® J b ® I c)® 

^\c{J r ® la® Jb® Ic) + ^‘bc{J r ® J a ® Ib ® Ic)® 
^ABciJr ®Ia®Ib® Ic) + ® la® Ib® Ic)- 

The covariance matrix when only A has fixed effects is: 

V{Y) = aR{Ir ®Ja®Jb® Jc) + {Ir ® la® Jb® Jc) + 
a\{Jr ® Ja® Ib® Jc) + ®Ja®Ib® Jc) + 

CrAB{Jr ®Ia®Ib® J c) + ^^Iab^Ic ®Ia®Ib® J c) + 

ac{Jr ®Ja®Jb® Ic) + (^Ac{Jr ® la® Jb® Ic) + 

^‘bc{J r ® j a ® Ib ® I c) + ^\bc{J r ® I a ® Ib ® Ic)® 

Crl^{Ir ®Ia®Ib® Ic)- 
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The covariance matrix when only B has fixed effects is: 


V (Y) = (T%{Ir ® Ja®Jb® Jc) + (y\{Jr J c) + 

CJg^ (-^r ® I a ® J b ® J ^cb ® J a ^ ^b ® J c)T 
(TABiJr ^la'^Ib'Si J c) + ‘S)Ia^Ib® J c) + 

ac{Jr ^Ja'SlJb'Sl Ic) + (^AciJr ^la'SlJb’^ Ic) + 

r ^ Ja 'S’ Ib ^ Ic) + ^\Bci'I r ® la'^ Ib® -^c) + 

®Ia®Ib® Ic)- 

The covariance matrix when only C has hxed effects is: 

V (V) = ajiilr 0 Ja^Jb^ Jc) + 0-\{Jr ®Ia®Jb® J c) + 

0-l^{Ir 0la'S)Jb'S> Jc) + 0-%{Jr ® Ja®Ib® J c) + 
CJg^ (yI-r®Ja®Ib®Jc) T ^ AB^J r ® I a ^ I b ® J c)~^ 
O-lAB^Ir ®Ia®Ib ®Jc) + f^AciJr ^Ia®Jb® Ic) + 

^“sciJ r ® Ja Ib ^ Ic) + ^“jiBciJ r ^ la Ib ^ -^c) + 
al^{Ir <^Ia®Ib<Si Ic)- 

The covariance matrix when only A has random effects is: 

V (Y) = a\{Ir ® Ja®Jb® Jc) + 0-\{Jr ®Ia®Jb® J c) + 

al^{Ir ® la ®Jb®Jc) + 0-lg{Ir 0Ja<S'Ib'^ J c) + 

(JABiJr <^Ia®Ib<S' Jc) + 0la<^Ib<S' Jc) + 

<^AciJ r ® la® Jb® I c) + ^“aBC^J r ® la ® Ib ® -^c) + 
al^ilr ®Ia®Ib® Ic)- 

The covariance matrix when only B has random effects is: 

V{Y) = (T\{Ir ®Ja®Jb® Jc) + <7g^ {Ir ® la® Jb® Jc) + 
aB{Jr ® Ja® Ib® Jc) + 0-lg{Ir ®Ja®Ib® Jc) + 

^ AB {J r ® I a ® Ib ® J c) Y ^cab ® ^ ^ ® ® Jc)Y 

^“sciJr ® Ja ® Ib ® Ic) + r ® la ® Ib ® -^c) + 

Crl^{Ir ®Ia®Ib® Ic)- 

The covariance matrix when only C has random effects is: 

V{Y) = aB{Ir ®Ja®Jb® Jc) + (^Ia {Ir ® la® Jb® Jc) + 
®Ja®Ib® Jc) + (^CAB^^r ®Ia®Ib® Jc) + 
ac{Jr ®Ja®Jb® Ic) + (^AciJr ® la® Jb® Ic) + 
^“bC^J r ® j a ® Ib ® Ic) + r ® la ® Ib ® -^c) + 

al^ilr ®Ia®Ib® Ic)- 
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3 Expected mean squares 


To illustrate how to obtain the expected mean squares E{MS) we will show the process for the 
random blocks R (For the remaining cases, since the procedure is similar, we will only present 
the final value without the respective development): First, take SSr in Tableland calculate its 
expected value: 


E{SSr) = abc ^ E{Rh - R. + - e^., + + cabh.. - eAB... + - ej )^, 

h=l 

since the product of errors and factors is always zero under expectation, E{SSb) equals 

r r 

abc'^E{Rh - R)"^ + abc'^E{eAu. -^a.. Egb^. -gb.. +eABh.. - gab... + 6 *^,,, -Gt....f, 

h=l h=l 

and since the errors are independent between themselves, 

r r r 

E{SSr) = abc^E{Rh -R.)^ + abc^EicA^ -eA..f + abcY,E{eB^ -eB..f 

h=l h=l h=l 


+ abc E E{gaBh.. - gab...)'^ + abc'^E{et^ - a... 


h=l 


h=l 


Therefore, taking into account that Ug = E{e^) — Fi^(e), and that E{e) =0 for every error in the 
model, 

^ _ (r - IW^ 

E{SSr) = abc'^EiRh - R.f + abc- -^ 


+ abc 


h=l 


+ abc- 


+ abc- 


b ab abc 

Then, taking E{SSr) and dividing it by its df, we obtain the expected mean square for hxed blocks: 

r 


abc 


E{MSr) = ^ ^ E{Rh -R.f + bcal^ + acal^ + 


2 I 2 


( 1 ) 


h=l 


Now, for the more interesting case of random blocks, we get: 

E{MSr) = abca\ + bcal^ + acaj^ + caj^^ + 

Note that E{MSr) will remain unchanged regardless the model we are considering. This is also 
true for the expectation of the mean square of each error involved. So we mention these here and 
will omit them in the particular description of the E{MSys for each model: 

E{MSeA) = ^Gal^ + cal^^ + a^, 

EiMSeg) = acalg + cal^^ + a^, 

E{MSe,,) = cal^^+al, 

E{MSe,) = al- 

Finally, note also that every interaction involving a random effect will be random. So in the 
following subsections, to avoid confusion, we present explicitly all the E{MSys for every model. 
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3.1 Expected mean squares for the fixed effects model 

When the effects are fixed (constant), by definition it is sufficient to suppress the expectation 
operator of the mean squares considered. Thus we get: 


JL. _ _ _ _ _ _ _ 

E(MSa) = - y^(Ai-A +ABi -AB + ACi - AC + ABCi -ABC f 

a — 1 

i=l 

+ bcal^+cal^^+al, 

h 

E{MSb) = Y^{Bj - - A^„f 

i=i 

+ acal^ + cal^^ + al, 


E{MSab) = 


cr 


a b 




^ '^{ABij - AB,j - ABi, + AB,, 


+ ABCij, — ABC ,j, — ABCi,, + ABC,,,) + ccTg^^ + Ug^, 
abr 


E{MSc) = ^ -C,+ AC,k - AC„ + BC,k - BC„ + ABC„k - ABC+ al, 


k=l 


E{MSac) = 


br 


- AC,k - ACi, + AC„ 


+ 'ABCi,k - 'ABC„k - ABCi,, + ABC,,,)^ + a, 

b c 


2 , ,2 
et ’ 


E{MSbc) = 


ar 


(6- l)(c- 1) 

^ > 7 = 1 k=l 




+ ABC,jk - ABC„k - ABC,j, + ABC„X + 

a b c 


E{MSabc) = 






jk 


- ABCi,k + ABC„k - ABCij, + ABC,j, + ABCi,, - ABC„y + cr, 


3.2 Expected mean squares for the random effects model 

When a factor, say A, has random effects, we will assume that the effects of A have distribution 
N (0, fj^). Then for the random effects model, the effects of A, B and C will be random, independent, 
and normally distributed with mean 0 and variance a\, a\ and respectively. The interactions 
AB, AC, BC and ABC will have normal distribution with mean 0 and variance 
and cr^BC^ respectively. We also assume that the effects are independent between them. So we get: 
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E{MSa) = bcra\ + bca^^ + cral^ + c(Tl^g,+bra\c + ra^BC + 
E{MSb) = acral + aca^^ + cral^ + cal^^,+aralc + raise + 
E{MSab) = era'll + ca^^^ + ral^fj + a^, 

E{MSC) = abral + bra\c + aral(^ + raise + *^et) 

E{MSac) = brale + raise + ^“It^ 

E{MSbc) = arale + raise + 

E{MSabc) = ral^fj + a^. 


3.3 Expected mean squares when only A is fixed 

In this case, B and C will be random with variances al and respectively. Also, AB, AC, BC 
and ABC will be random with variances al^, al(j, al^j and al^fj, respectively. So we obtain: 


E{MSa) = —— 22{Ai - Af + bcal^ + cral^ + cal^^ + brale + ral^e + 

® ^ *=i 

E{MSb) = acral + aca^^ + crals + + arale + raise + 

E{MSab) = crals + ca^^^ + raise + cr^, 

E{MSe) = abral + brale + arale + raise + '^et> 

E{MSac) = bralu + ral^Q + a^, 

E{MSbc) = arale + raise + '^eti 
E{MSabc) = ral^e + '^et- 


3.4 Expected mean squares when only B is fixed 

Here we have that A and C will be random with variances al and cr^, respectively. Also, AB, 
AC, BC and ABC will be random with variances al^, al^j, al^j and al^^, respectively. So we 
obtain: 


E{MSa) = hcral + bca^^ + crals + ca^^g + brale + raise + 

b 

E{MSb) = ^ ^ + ®ccTg^ + crals + ca^^^ + arale + raise + 

i=i 

E{MSab) = crals + ca^^^ + raise + cr^, 

E{MSe) = abral + brahj + aral(j + ral^e + 

E{MSAe) = brale + calse + c^l, 

E{MSbc) = arale + raise + '^et) 

E{MSabc) = raise + <^et' 


7 




3.5 Expected mean squares when only C is fixed 

In this case, A and B are random with variances <7^ and respectively. Also, AB, AC, BC and 
ABC will be random with variances ^\bc^ respectively. So we obtain: 


E{MSa) = hcra\ + hcal^ + cra\s + + bra\c + rcrABC + 

E{MSb) = acra^ + + cra\B + cal^^ + ara%Q + ra\BC + 

E{MSab) = CraAB + ''^^abc + 

E{MSc) = ^ ^ + '^^\bC + ^“eti 

^ k=l 

E{MSac) = bra\c + f^ABC + 

E{MSbc) = CL^Cr'^C + ’’’^ABC + '^et) 

E{MSabc) = ^(^ABC + <^64 • 

3.6 Expected mean squares when only A is random 

In this case, A, AB, AC and ABC are random with variance a\, cr\^, a^Q and (t\bc^ respectively. 
Thus, we get: 


E{M Sa) = hcra\ + bcal^ + cra\B + + bra\c + rcr^BC + 

b 


E{MSb) — ^ ^ — B, + BCj, — BC,)^ + accjg^ + cra\B + 

i=i 

E{MSab) = cra\s + ca^^g + ra^sc + 

E{AISc) = -—- ~ C, + BC k — BC + bra\(j + raABC + 


+ '>'(^ABC + '^et; 


E{MSac) — bra AC + '^^\bc + 

b c 

E{MSbc) = E Y.(BC,k - ^.k - + ra\BC + 

^ ’ i=i k=i 

E{MSaBc) = I^CTaBC + ^et- 


3.7 Expected mean squares when only B is random 

In this case, B, AB, BC and ABC are random with variance a‘^, a\^, a\(j and cr^ec’ respectively. 
Thus, we get: 
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E{MSa) 

E{MSb) 

E{MSab) 

E{MSc) 

E{MSac) 

E{MSbc) 

E{MSabc) 


“ _ _ _ 

_ I + ACi, — + &C(Tg^ + craAB + '’"^abc + ^et) 

1 = 1 

acral + aca^^ + cral^ + c^^Iab + 
crCTAB + (^^BAB + ’’"^ABC + ^Bt ) 

^ _ _ _ 

_ ^ — C, + AC,k — AC,)"^ + aralfj + raABC + *^64> 

fe=i 

, a c 

7—_ ^ y^lACjk — AC A — ACi, + AC,,)'^ + raABC + *^ 64 ) 

4=1 k=l 

+ ^^\bC + <^641 
’’’^ABC + <^64 • 


3.8 Expected mean squares when only C is random 

Here C, AC, BC and ABC are random with variance a^, cr\c'’ ^BC ^ABC^ respectively. Thus, 
we get: 


“ _ _ _ 

E{MSa) = ^ ^ — A, + ABi, — AB,)'^ + boa'll + co-g^^ + hra\c + i^^abc + <^64) 

i=l 

b 

E{MSb) = ^ ^ '^^(Ej — B, + AB j — AB )‘^ + acal^ + cUg^^ + aral(j + ra^BC + *^64) 

j=i 

a b 


E{MSab) — (-Q _ _ 1) 'Yl - AB,j - ABi, + AB,)"^ + ra\BC + ^^Iab + '^e 4 > 

E{MSc) = abral + bralc + ^.ralc + '^'^ABC + '^641 

E{MSac) = bra\c + ’’’^Ibc + *^ 64 ) 

E{MSbc) = + '''^ABC + *^64) 

E{MSabc) = f^lBC + '^64- 


4 F tests 


In this section we present the T tests in tables. When required, we wil l specify the approximated 


df by means of the famous estimator developed by Satterthwaite ( 1946l i: for the cases in which the 


co mplex estimation is a funct ion of two variance components we will use the estimator proposed 
bv lAmes and Webster (119911 1. which is a correction to Satterthwaite for this particular case. Let 
us start with the estimator by Satterthwaite: 

If 9 is variance which is a linear combination of m independent variances, i.e., if 0 = o-iOi, 
with estimator 6 = YFLi we say that 0 is a complex estimator of 9. Since for our case the 

coefficients ai = 1, for i = 1,..., m, we will omit them on what follows. For the cases in which the 
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variance estimator is complex, Satterthwaite (194^) proposed the following estimator for the df: 


fs = 


{TZi 

ET=iMSf/n.^ 


( 2 ) 


where rii are the df of the source of variation corresponding to i. This is so because ^ can be 
approximated to a with / degrees of freedom. This estimato r is widely known, so we proceed 
to present the estimator proposed bv lAmes and Webster (|l991h : 


When the variance 0 is a function of two variance components 9i and 62 -, call (pi = 1 and 
4>2 = 92 IO 1 , and consider the class of estimators given by <1)2 = rMS' 2 /MS'i, where r is a constant, 
then we can approximate the df by 


fawir) = 


(ELi - 


Et=i 


(3) 


Note that faw{^) = fs and that min(ni,n 2 ) < fawii") < ni + n 2 - Thus, we can vary r in order to 
get better properties. For instance, 

^ 712 f 2(ni +772-2) 

T = - I -h 1 

772 - 2 V ni(r72 - 4) 

minimizes the mean square of the error of l/(j)2- Also r* > 1 and fawii"*) < fs- In this paper, every 
time we calculate the Ames-Webster estimator ([3]), we will also calculate its respective value r*. 
Using the Ames-Webster approach we have two possible estimations for every value of r. Then, if 
both of them are less than fg, it is advisable to use the larger one, since the smaller one usually 
has a negative bias. 

With these tools at hand, we proceed to present the F test for every model. The hrst column in 
each of the following tables will be the source of variation, the second one will tell us if the effects 
are random or fixed, the third one will be the corresponding F test and the last one will be the null 
hypothesis under consideration. When the effects are random, the null hypothesis will be that the 
corresponding variance of the source has 0 variance; when the effects are hxed, the null hypothesis 
will be that all effects are equal (to 0). 


4.1 F tests vi^hen all effects are fixed 


To construct the F tests in Table [21 we use the expected mean squares found in Subsection 13.11 
Note that R will have the same structure for the F test, regardless of it being constant or random 
(although, of course, the hypothesis will change). 

Using the Satterthwaite estimator in ([2]), we approximate the df for R as: 


Vl 


V 2 


{MSr + MS, 


SABJ 


MS% 
r—1 


+ 




(r-l)(a-l)(6-l) 


MSi^ MSl^ ’ 

+-1)(«-1) + +-!)+-!) 


where vi and V2 are the df for the enumerator and the denominator, respectively. 

When we adjust using the Ames-Webster estimator (|3]), we obtain two estimators for each case. 
First let us see the df for the numerator: Let MSi = MSr and MS2 = MS^^b^ then 


10 






















Table 2: F tests for the fixed effects model 


Source 

Effect 

F 

Ho 

R 

/ 

M bR-\-M tie. ^ 
MSe^+MSer, 

0 

II 

b 

A 

/ 

M tb'A 

MSeA 

II 

to 

II 

II 

s 

II 

0 

eA 

r 

MSe^ 

MS. 

0 

II 

b 

B 

/ 

MSb 

MS.„ 

Bi = B2 = • ■ • = Bb = 0 

e-B 

r 

MS.„ 

<=o 

AB 

f 

MSab 

MS.ab 

{AB),j = 0 , Vi, Vj. 

eAB 

r 

MS.^b 

MSet 

= 0 

C 

f 

MS., 

II 

II 

II 

P 

II 

0 

AC 

f 

M bAC 

MS., 

{AC),k = 0 , Vi, Vfc. 

BC 

f 

MbBC 

MS., 

{BC),k = 0 , Vj, Vfc. 

ABC 

f 

Mb Ann 

MS., 

{ABC)„k = 0 , Vz, Vj, Vfc. 

et 

r 

- 



(r-l)(a-l)(b-l) / 2[{r-l){a-l){b-l) + r-3] A 

(r - l)(a - 1)(5 - 1) - 2 V(r - l)[(r - l)(a - l){b - 1) - 4] ^ ^ ’ 


~ , ^_{l+p,MSe^jMSRf 

Jaw[Pl) ^ (piM5e^„/MSa)^ ’ 
r—1 ' (r—l)(a—l)(b—1) 


on the other hand, when MSi = MSe^^ and MS 2 = MSr: 

* ^ r-l f 2[{r - l)(a - l)(b - 1) + r - 3] 
r —(r — l)(a — 1)(6 — l)(r — 5) 



1 _ (plMSH/MSe^^f ■ 

(r—l)(a—l)(b—1) ' r—1 

Now, for the denominator, when MSi = MSe^ and MS 2 = MSeg we have: 

^ (r-l)(b-l) / 2{(r-l)[(a-l) + (b-l)]-2} \ 

(r-l)(6-l)-2 V(r-l)(a-l)[(r-l)(6-l)-4]^ )' 


faw{P2) = 


il+P2MSejMS, 


SAJ 


1 


(r-l)(a-l) 


+ 


{p2MSe„/M-3)^ ’ 
(r-l)(b-l) 


and when MSi = MSeg and MS 2 = MSsa we have: 


P2 = 


(r - l)(a - 1) / 2{(r - l)[(a - 1) + (6 - 1)] - 2} 


(r — l)(a — 1) — 2 \{r — 1)(6 — l)[(r — l)(a — 1) — 4] 


+ 1 


faw{P 2 ) = 


{l+p*^MSe^/MS, 


es) 


(r—1){6— 1 ) (r—l)(a—1) 

The estimators for the df of R will always be the same. For this reason they will be omitted on 
the tables to come. 


(p*M5e./M-5)2 
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4.2 F tests when all effects are random 

When all effects are random, we construct the F tests on Table [3] based on the mean squares 
developed in Subsection 13.21 


Table 3: F tests for the random effects model 


Source 

Effect 

F 

Ho 

R 

r 

AlitR + MBej^g 

MSe^+MSe„ 

4 = 0 

A 

r 

M S>a + M -\-MtiABC 

MSe^+MSAB+MSAC 

II 

o 

eA 

r 

M 

o 

II 

b 

B 

r 

+Ai6’ASC 

MSej:t-\-MSAB-\-MSBC 

4 = 0 

e_B 

r 

MSes 

<^BR = 0 

AB 

r 

^AB 

MSe 4 R -\-MSabc 

^AB = 0 

vab 

r 

MSe, 

= 0 

C 

r 

A1 !Sc-\-Af k^abc 

MS Ac+MSnn 

II 

o 

AC 

r 

aISac 

AISabc 

o 

II 

o 

b 

BC 

r 

M SbC 

AISabc 

^BC ~ 0 

ABC 

r 

MS ABC 

AlSe, 

’^ABC — 0 

et 

r 

- 



Since the complex estimators for effects A and B in Table [3] have three variance components, 
we will use only ([2]) with them to find their approximate df. For the effects of A, the df in the 
numerator and the denominator vi and V 2 , respectively, will be given by: 


Vl 


V2 


{MS A + MSe^^ + MS ABC? 

MS\ I 

a—1 ' (a—1)(6—l)(r—1) ' (a—1)(6—l)(c—1) 

{MSe^ + MSab + MSAcf 

MS\g MS\„ 

(r-l)(a-l) (a-l)(fe-l) + (a-l)(c-l) 


(4) 


With B, the df of the F test will be for the numerator and denominator respectively: 


Vl 


V2 


{MSb + MSe^^ + MSabc? 

MSI I I 

-b^ + (a-l)(6-l)(r-l) + (a-l)(fe-l)(c-l) 


{MSe^ + MSab + MSbc? 


(r-l)(b-l) 


+ 


(a-l)(6-l) 


+ 


MS%a 

(6-l)(c-l) 


(5) 


The df for AB approximated by ([2]) will be respectively for the numerator and the denominator: 


{MSab + MSe,f 

I ’ 

(a—1)(6—1) ' a6(c—l)(r—1) 

{MSe^^ + MS ABC? 

(a-l)(b-l)(r-l) + (a-l)(6-l)(c-l) 


( 6 ) 
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Adjusting by means of ([3]), there are two possible estimators in each case. First, let us see the 
degrees of freedom in the numerator. Let MSi = MSab and MS2 = MSe^., then 


Pi 


ab{c — l)(r — 1) / 2[(a 

ab{c — l)(r — 1) — 2 \ (o 


1)(6 — 1) + ab{c — l)(r — 1) — 2] 
1)(6 — l)[a6(c — l)(r — 1) — 4] 


faw {Pl) 


{l+PlMSejMSAB? 

1 , {piMS^JMSabY ’ 

(a—1)(6—1) ' a6(c—l)(r—1) 


on the other hand, when MSi = MSe^ and MS 2 = MSab'- 


, (a-l)(6-l) /2[(a-1)(5-1)+a6(c-l)(r-1)-2] 

“ (o - 1)(6 - 1) - 2 V ab{c - l)(r - l)[(a - 1)(6 - 1) - 4] 


(7) 


faw{p\) -- 

And for the denominator, when MSi 


{l+plMSAB/MSe,? 

1 , (p*MSAB/MSet)‘^ ■ 

afo(c—l)(r—1) ' (a—1)(&—1) 

: MSej^g and MS2 = MSabc we have: 


( 8 ) 


(g- 1)(6- l)(c- 1) / 2[(a-l)(6-l)(c + r-2)-2] 

(g - 1)(6 - l)(c - 1) - 2 V (o - 1)(6 - l)(r - l)[(g - 1)(6 - l)(c - 1) - 4] 


faw{P2) = 


{I+P 2 MSABC/MS, 


^AB) 


+ 


(p2MSABc/MSej^g)^ ’ 


(a-l)(b-l)(r-l) (a-l)(b-l)(c-l) 

when MSi = MSabc and MS2 = MSe^g we get: 

(g-l)(6-l)(r-l) ( 2[(g - 1)(6-l)(c + r - 2) - 2] 


P 2 = 


(g - 1)(6 - l)(r - 1) - 2 V (a - l){b - l)(c - l)[(g - l){b - l)(r - 1) - 4] 


+ 1 


(9) 


faw{P*2) = 


(1 +p1,MSe^g/MSABC? 


+ 


(p-^MSe^g/MSABc)^ 


With C, using /*, the df for the numerator and denominator will be respectively: 

{MSc + MSabc)"^ 


Vl = 


V2 = 


MSI I MSlg^ ’ 
c—1 ' (a—l)(fe—l)(c—1) 

{MSac + MSbc)'^ 


(a-l)(c-l) ^ (a-l)(6-l)(c-l) 


( 10 ) 


( 11 ) 


For faw these were the estimators for the df of the numerator when MSi = MSc and MS 2 = 
MSabc- 


(g- 1)(6- l)(c- 1) / 2[(g-l)(6-l)(c-l) + c-3] \ 

(g - 1)(5 - l)(c - 1) - 2 V(c - l)[(g - 1)(6 - l)(c - 1) - 4] + 
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~ {I+PiMSabc/MSc? _ 

jypi) 1 , [piMSabc/MScY ’ 

c—1 ‘ (a—l)(b—l)(c—1) 

still for the numerator, but exchanging the order of MSi and MS'2? we obtain: 


* ^ c- 1 ( 2[(a-l)(6-l)(c-l) + c-3] \ 

c —3\ (o — 1)(6 — l)(c — l)(c — 5) /’ 


(12) 


m) 


(1 + pIMSc/MSabc)'^ 

1 , {p*,MSc/MSABcr • 

(a—l)(fe—l)(c—1) c—1 


(13) 


For the denominator, taking MSi = MSac and MS2 = MSbc-, we get the following estimations: 


(6-l)(c-l) / 2{(c-l)[(a-l) + (b-l)]-2} \ 

(6 _ l)(c - 1) - 2 V (a - l)(c - 1)[(6 - l)(c - 1) - 4] ^ ) 


f{P2) 


{1 + P 2 M S BC / M s Ac)^ 

1 I {P2MSbC/^^Sac)'^ ’ 

(a-l)(c-l) + (6-l)(c-l) 


(14) 


once again for the denominator, but exchanging to MS\ = MSbc and MS 2 = MSac, we get: 

(a-l)(c-l) / 2{(c-l)[(a-l) + (fe-l)]-2} \ 

P2 («_i)(c_i)_2 V(6-l)(c-l)[(a-l)(c-l)-4]+ b 



{1+P*2MSac/MSbc? 

1 , {p%mSac/mSbcY • 

(h-l){c-l) + (a-l)(c-l) 


(15) 


4.3 F tests when only one factor has fixed effects 

With respect to Table [3l the only difference for the three cases considered here (only A has fixed 
effects, only B has fixed effects, and only C has fixed effects) will occur in the row corresponding 
to the fixed effect: first, obviously, its effect will be / instead of r; second, its null hypothesis will 
be about the equality of all treatments inside that factor. So when A is the only factor of fixed 
effects, its effect is / and its null hypothesis is Ai = • • • = ^^ = 0, all other fields remaining equal 
to Table [Sj when B is the only factor with fixed effects, its effect is / and its null hypothesis is 

= • • • = = 0, all other fields remaining equal to Tabled and when the only fixed effects are 

those corresponding to C, its value at effect is / and the null hypothesis will be Ci = • • • = Cc = 0, 
all other fields remaining equal to Table [3l This can be easily verified with the information in 
Subsections 13.3113.41 and 13.51 

Since, in particular, the structure of the complex variance estimators is identical to the structure 
of the model with random effects, the approximate df for each of these three cases are exactly the 
same to those found in Subsection 14.21 
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Table 4: F tests when only A is random 


Fuente 

Efecto 

F 

Ho 

R 

r 

IVI Bg + lVI 
MSe^+MSe„ 

0 

II 

b 

A 

r 

IVl iV'i ib ABC 

MS^^+MSab+MSac 

0 

II 

b 

eA 

r 

MScab 

0 

II 

b 

B 

f 

MS^b+MSab 

Bi = B 2 = ■ ■ ■ = Bb = 0 

e_B 

r 

MS.ab 

<=o 

AB 

r 

MbAB+M!^et 

MSe^n +^Sabc 

<^AB = 0 

eAB 

r 

MSTFe 

MSe, 

= 0 

C 

f 

i\l be 

MS AC 

II 

II 

II 

P 

II 

0 

AC 

r 

M bAC 

MS ABr 

^AC — 0 

BC 

f 

M bBC 

MS ARC 

(BC),fe = 0,Vj,fc. 

ABC 

r 

M bABC 

MSe, 

^ABC = 0 

e-t 

r 

- 



4.4 F tests when only A has random effects 

When the effects of A are random, we obtain Table 0] based on the -E(MS')’s found in Subsection 

ITGI 

The estimators for the df of A are those in Q. The estimators of the df by Satterthwaite for 
AB are those in Q; the estimators by Ames-Webster are given in equations d?]) and ([8|) for the 
numerator, and ([9]) and (llOjl for the denominator. Now we precede to evaluate the df for the F 
test of B, first by means of the Satterthwaite estimator in equation ([2]): 


vi 


V 2 


{MSb + MS, 


SAB) 


MSI I 

b—1 ' (a—l)(fo—l)(r—1) 

(MSe^ + MSab? 


MSF 


MSls 


_ ^B I ^AB 

(6-l)h-l) (a-l)(6-l) 


(16) 


Still with B, the first Ames-Webster estimator for the df of the numerator of the F, taking MSi = 
AISb and MS 2 = MSe^g will be: 


(a-l)(6-l)(r-l) /2 {{b-l)[{a-l){r-l) + l]-2} \ 

(a_l)(5-l)(r-l)-2 V(&-l)[(a-l)(6-l)(r-l)-4] + J’ 


~ (l+p^MSe^g/MSB? 

FPV ^ ^ (piMS.^g/MSgy ’ 

b-l (a-l)(bTi)(r-l) 

and exchanging the order to MSi = MSe^^g and MS 2 = MSb, we obtain: 

* ^ / 2{(6-l)[(a-l)(r-l) + l]-2} \ 

6-3 V (a-1)(5-l)(r-l)(6-5) ) 


f{p\) 


{l+plMSB/MSe^gf 
1 , {plMSg/MS.^^f 

(a—l){'r—1) ' 6—1 
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For the denominator of the the F test of B, taking MSi = MSg^ y MS 2 = MSab'- 

(a-l)(6-l) / 2[(6-l)(a + r-2)-2] \ 

(a-l)(6-l)-2 V(6-l)(r-l)[(a-l)(6-l)-4]+ )' 


f(.P2) 


{l+P2 MSAB/MSe^)^ 

1 {p2MSABlMS,g)^ ’ 

{b-l){r-l) + (a-l)(6-l) 


and exchanging the order of MSi and MS' 2 : 


(b-l)(r-l) / 2[(6-l)(Q + r-2)-2] 

P2 (5 _ i)(, _ 1) _ 2 V (a - 1)(6 - 1)[(6 - l)(r - 1) - 4] 


f{P*2) 


{l+p*2MSeg/MSAB)^ 

1 , {plMS.glMSAB? ■ 


4.5 F tests when only B has random effects 

When only A and C have fixed effects, based on Subsection (|3.7I) . we get the following: 


Table 5: F tests when only B is random 


Source 

Effect 

F 

Ho 

R 

r 

MSe^+MSe„ 

4 = 0 

A 

f 

iVL 

MSe^+MSAB 

II 

to 

II 

II 

II 

0 

eA 

r 

MSbab 


B 

r 

MiSB+M!5eji^^Q +M1 :Sabc 
MS e„+MSAB+MSBC 

0 

il 

C4 CQ 

b 

e-B 

r 

MSeg 

MSbab 

= 0 

AB 

r 

-^‘5'e 4 R +MSaBC 

’^AB = 0 

eAB 

r 

MS., 

= 0 

C 

f 

M be 

MSnn 

II 

II 

II 

P 

II 

0 

AC 

f 

M Sac 

MSabci 

{AC)ik = 0, Vz, k. 

BC 

r 

M bBC 

MS ABC 

^BC = 0 

ABC 

r 

Mb ABC 

MS., 

^ABC = 0 

et 

r 

- 



The approximated df for the F test of B were found using the Satterthwaite estimator (j5]) . 
The approximation of the df for AB using Satterthwaite is given by ([6]); using Ames-Webster, the 
estimator for the df of AB are given in equations d?]) and dH) for the numerator, and dH) and (uni) 
for the denominator. Now we precede to evaluate the df for the F test of A, first by means of the 
Satterthwaite estimator in equation dl]): 


vi 


V 2 


{MSA + MSg^gf 

MSI I ^^^B 
a—1 (a—1)(6—l)(r—1) 

{MSg^ + MSab? 

I Mslg 
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The Ames-Webster estimator for the numerator is the following when MSi = MSa and MS 2 = 


(a-l)(b-l)(r-l) / 2 {{a-l)[{b-l){r-l) + l]- 2 } \ 

(a - 1)(6 - l)(r - 1) - 2 V (« - l)[(a - 1)(6 - l)(r - 1) - 4] + ^ ’ 


, ^ { 1 +PiMS,^JMSa)\ 

Jm) 1 , (piMS,^„/MSa)^ ’ 

a—1 (a—1)(6—l)(r—1) 


still with the numerator but taking MSi = MSe^g and MS2 = MSa, we get: 

^ ^ / 2{(a-l)[(b-l)(r-l) + l]-2} \ 

^ a —3\ (a — 1)(5 — l)(r — l)(a — 5) / 


fipl) 


(l+plMSA/MSe^g)^ 

1 , (plMSA/MSe^g)^ ■ 

(a—l)(fe—l)(r—1) w fl,_i 


For the denominator, doing MSi = MSe^ and MS2 = MSab, we get: 


(g- 1)(6- 1) / 2[(a-l)(6 + r-2)-2] 

(g - 1)(6 - 1) - 2 V(g - l)(r - l)[(g - 1)(6 - 1) - 4] 


f{P2) 


(1+P2M5AS/M5e^)2 
1 , (p 2 MSabIMS,^Y ’ 

(a-l)(r-l) w (a-l)(b-l) 


finally, exchanging the order of MSi and MS 2 , we obtain: 


(g- l)(r- 1) / 2[(g-l)(fe + r-2)-2] 

(g - l)(r - 1) - 2 V(g - 1)(6 - l)[(g - l)(r - 1) - 4] 


{l+P* 2 MS,^/MSABf 

^ ^ (plMSe^/MSAB)^ ' 

(a-m-l) + (a-l)(r-l) 

4.6 F tests when only C has random effects 

Table [6] below was constructed using the E{MSys in Subsection 13.81 Note that the approximate 
degrees of freedom for AB were described in equations ([6j) by means of Satterthwaite; also for AB, 
the approximations of its degrees of freedom using Ames-Webster were given in equations ([7]) and 
(|8|) for the numerator , and Q and m for the denominator. 

For C, its approximate df using Satterthwaite were found in (jllj) . And the Ames-Webster 
estimators of the df of C are given by (I12p and ()13p for the numerator, and by ()14p and (IlSp for 
the denominator. 

We precede to evaluate the approximate df for the F test of A, hrst by means of the Satterth¬ 
waite estimator in equation ([2]): 
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Table 6: F tests when only C is random 


Source 

Effect 

F 

Ho 

R 

r 

M bR-\-M *^€^3 

-^1‘S'e ^ -\-M.Se. R 

0 

II 

b 

A 

f 

iW h 

MSeji +MSaC 

II 

II 

II 

0 

II 

0 

ga 

r 

MSe 4 f, 

<=o 

B 

f 

M 6's+-/Vi i>'et 

MSe„+MSBC 

Bi = B 2 = ■ • • = Bb = 0 

e-B 

r 

JUS'ej, 

MSbab 

< = 0 

AB 

f 

1\1 + Ai ib'ey. 

-^‘5'e 4 R +MSaBC 

{AB)ij = 0, Vi,j. 

gab 

r 

MS,, 

^e.B=0 

C 

r 

M ^ABC 

MSac.+MSbc 

0 

II 

b 

AC 

r 

MSac 

MSa Bn 

0 

II 

0 

b 

BC 

r 

M Sbc 

MS A Bn 

^BC — 0 

ABC 

r 

Mb ABn 

MS,, 

^ABC = 0 

et 

r 

- 



(MSA + MSe,)^ 

^ MSI ’ 

a—1 ' ab{c—l){r—l) 

{MSe^ + MS AC? 

MSj^ , MS^r: 

(r-l)(a-l) (a-l)(c-l) 


The Ames-Webster estimator for the numerator is the following when MSi = MSa and MS2 

MSe,: 


Pi 


ab{c — l)(r — 1) 
ab{c — l)(r — 1) — 2 


( 2[ab{c - l){r - 1) + a - 3] \ 


~ {l+piMSejMSA)\ 

J[Pl) 1 (piMSeJMSA)^ ’ 

a—1 ' a6(c—l)(r—1) 

still with the numerator but taking MSi = MSet and MS2 = MSa, "we get: 


Pi = 


a — 1 
a — 3 


/ 2[a6(c — l)(r — 1) -b o — 3] 
\ a6(c — l)(r — l)(a — 5) 


m) 


{l+plMSA/MSe,? 

1 , {p*MSa/MS,,F ' 

ab{c—l){r—l) a—1 


For the denominator, doing MSi = MS^^ and MS 2 = MSac, we get: 


(a-l)(c-l) / 2[(o-l)(r + c-2)-2] 

(« _ i)(c _ 1) _ 2 V(a - l)(r - !)[(« - l)(c - 1) - 4] 
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f{P 2 ) = 


{l+P 2 MSAc/MSe^)^ 

1 , {p2MSAclMSe^f ’ 

(a-l)(r-l) + (a-l)(c-l) 


finally, exchanging the order of MSi and MS 2 , we obtain: 


P2 = 


(a-l)(r-l) / 2[(a-l)(r + c-2)-2] 

(a - l)(r - 1) - 2 V(a - l)(c - l)[(a - l)(r - 1) - 4] 

(1 +p*^MSejMSAC? 


+ 1 


f{p*2) = 


1 , {pIms.^/mSacY' 

{a-l)[c-l) + (a-l)(r-l) 


Now, we evaluate the approximate df for the F test of B, first by means of the Satterthwaite 
estimator in equation Q: 


Vl = 


V 2 = 


{MSb + MSe,f 


MSI, 


_ 

6—1 ' ab{c—l)(r—l) 

{MSe^ + MSbc? 

MSi^ MSl^ 

(r-l){b-l) + (6-l)(c-l) 


The Ames-Webster estimator for the numerator is the following when MSi = MSb and MS 2 = 
MS,- 


Pi = 


ab{c—l){r — l) f 2[a6(c — l)(r — 1)-|-6 — 3] 


ab{c — l)(r — 1) — 2 V(^ “ l)[a6(c — l)(r — 1) — 4] 

{ 1 +PlMSejMSB? 


+ 1 


fiPi) = 


+ 


(piMSeJMSgY 


b—1 ' ab{c—l)(r—l) 

still with the numerator but taking MSi = MSg^ and MS 2 = MSb, we get: 

b — 1 f 2[a6(c — l)(r — 1) -b 6 — 3] 


Pi = 


6-3 

fipl) = 


ab{c — l)(r — 1)(6 — 5) 

{l+plMSB/MSe,f 

(plMSe/MSe-tY 

6-1 


+ 1 


_1_ I 

a6(c—l)(r—1) ' 

For the denominator, doing MSi = MSgg and MS 2 = MSbc, we get: 


P2 = 


(6-l)(c-l) 


2[(6-l)(r + c-2)-2] 


(6 _ l)(c - 1) - 2 V (6 - l)(r - 1)[(6 - l)(c - 1) - 4] 
f{P2) = 


+ 1 


{l+P 2 MSBc/MSe 




1 , {p2MSbcIMS,^Y ’ 

{b-l)(r-l) + (fe-l)(c-l) 

finally, exchanging the order of MSi and MS 2 , we obtain: 


* 

P2 = 


(6-l)(r-l) 


2[(6- l)(r + c-2) -2] 


(6 - l)(r - 1) - 2 V (6 - l)(c - 1)[(6 - l)(r - 1) - 4] 
{l+p^MSejMSBC? 


+ 1 


f{pl) = 


1 , (pIMS^JMSbcY 

(6-l)(c-l) (6-l)(r-l) 


19 



































5 Application 


In Izimmermanni (2004), a real example was considered when all the effects are fixed. The data in 
Tables [7] and [8] show the weight of 100 beans obtained by Luis Fernando Stone and Regis Vilela 
Bagatini on an experiment in 1998. It is a completely randomized design with two blocks on which 
each horizontal strip corresponds to the water layer irrigated, the vertical strips are soil tillage 
systems and the subplots are Nitrogen doses. The experiment was done at the Capivara farm in 
Embrapa Rice and Bean. 



The water layers (the vertical strips, A) are averaged irrigation levels as follows: 366.1 mm for 
the first horizontal strip, 335.1 mm for the second one, 315.7 mm for the third one, and 293.7 mm 
for the last one. There are three ways to prepare the soil (the vertical strips B)\ heavy harrowing 
for the hrst vertical strip, moldboard plowing for the second one, and notillage on the last one. The 
Nitrogen subdoses (C) inside the subplots are, respectively for each subplot, 0, 20 and 40 kg ha“^. 

The MS and the df needed to construct the F tests are shown on Table [9l These results show 
that there are significant effects on the water layer, its interaction with the soil, and the interaction 
of the three factors. 


Table 9: Examp 


e 


Source 

df 

MS 

F 

R 

1 

9.4758 


A 

3 

10.9903 

26.04 

SA 

3 

0.4220 


B 

2 

7.3937 

2.91 

es 

2 

2.5387 


AB 

6 

11.2718 

35.89 

e-AB 

6 

0.3141 


C 

2 

3.1476 

2.11 

AC 

6 

2.3759 

1.59 

BC 

4 

1.8678 

1.25 

ABC 

12 

3.2911 

2.21 

et 

24 

1.4921 



Now, in principle, it would be possible to analyze these data using other factorial experiments. 


20 









































Let us consider two of them: a three-way factorial experiment and split-split plot design. For these 
experiments, we present respectively Tables [TOl and [IT] below with the analysis of variance. 

The first obvious thing is that the B source (the soil) is not signihcative under the strip-split plot 
design, but is is signihcative under the two new models. It is important to note also the differences 
on the degrees of freedom for each model. In the strip-split plot design, for instance, the error of 
the effect B is represented by the interaction of the blocks and the soil, with 2 degrees of freedom; 
in the case of the factorial design the experimental error has 35 degrees of freedom (assuming the 
hxed effects model); and in the split-split plot design the error has 8 degrees of freedom, since the 
interaction blocks-layer-soil also includes the blocks-soil interaction. 

Now, these two cases are considered only for illustration purposes. However, it seems appro¬ 
priate to remember that the design of the held experiment dehnes the mathematical model and, 
therefore, the ANOVA. The most common mistake in statistics analyses applied to agriculture is 
to disrespect this design. This simple exercise shows that it is possible to carry on several analyses 
for a given data set, but only one is correct —the one following the held scheme. This held scheme 
is also determining the way to execute the randomization of treatments. 


Table 10: Factorial design 


Source 

df 

MS 

F 

R 

1 

9.4758 


A 

3 

10.9903 

8.73 

B 

2 

7.3936 

5.88 

AB 

6 

11.2718 

8.96 

C 

2 

3.1476 

2.50 

AC 

6 

2.3759 

1.89 

BC 

4 

1.8678 

1.48 

ABC 

12 

3.2911 

2.62 

et 

35 

1.2582 



Table 11: Split-split plot design 


Source 

df 

MS 

F 

R 

1 

9.4758 


A 

3 

10.9903 

7.37 

CA 

3 

0.4220 


B 

2 

7.3937 

4.96 

AB 

6 

11.2718 

7.55 

CAB 

8 

0.8702 


C 

2 

3.1476 

2.11 

AC 

6 

2.3759 

1.59 

BC 

4 

1.8678 

1.25 

ABC 

12 

3.2911 

2.21 

et 

24 

1.4921 
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